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ABSTRACT 

We consider a scenario where supermassive black holes form through direct accumulation 
of gas at the centre of proto-galaxies. In the first stage, the accumulated gas forms a super- 
massive star whose core collapses when the nuclear fuel is exhausted, forming a black hole 
of A/bh ~ IOOA/q. As the black hole starts accreting, it inflates the surrounding dense gas 
into an almost hydrostatic self-gravitating envelope, with at least 10 — 100 times the mass 
of the hole. We find that these "quasistars" suffer extremely high rates of mass loss through 
winds from their envelopes, in analogy to very massive stars such as 77-Carinae. Only for 
envelope masses greater than 2.8 x 10^(AfBH/100MQ)^/^^ is the envelope evaporation time- 
scale longer than the accretion time-scale of the black hole. This relation thus constitutes a 
"threshold growth line" above which quasistars can grow their internal black holes. Accretion 
rates can be 10 to 100 times the Eddington rate. The quasistars born in this "growth region" 
with 10^ — 10* A/q can grow black holes with masses between 10** to 10^ A/q, before crossing 
the threshold growth line and dispersing their envelopes in less than 10^ yr. This scenario 
therefore predicts that massive black hole seeds can be found only in dark matter halos with 
total masses larger than about IO^A/q, which can provide sufficiently high accretion rates to 
form such massive quasistars. 
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1 INTRODUCTION 

The formation of supermassive black holes (SMBHs) requires the 
inflow of > l(f Mq of gas from galactic scales into a region much 
smaller than a par sec. The presen ce of bright (~ 10''^ erg s~^) 
quasars at 2; > 6 JFan et aLluOOlh suggests that, at least in those 
cases, this inflow must have started at an earlier epoch and pro- 
ceeded at sufficiently high rate to have allowed the assemblage of 
black holes (BHs) of 10''i\fQ, in less than a Gyr. In proto-galactic 
dark matter halos with mass > 10^A/q, gas can be funnelled to- 
wards the centre at a rate of > Mq yr~^. 

If substantial frag mentation of gas into stars can be avoided 
(e.g.. lWise et alj|2008h . the collection of gas at such high rates can 
lead to the formation of a light (< lOOAfQ) black hole (BH) via 
direct collapse of gas or through the intermediate stage a supermas- 
sive star, whose core collapses to form a BH (Begelman 20\O^. The 
feedback from this "embryo" BH becomes important and it mod- 
ifies the structure of the pregalactic discs in its innermost « 100 
AU. The BH luminosity inflates the highly opaque gas in its vicin- 
ity into a pressure supported envelope, at least 10 — 100 times more 
massive than the BH itself. While the envelope is being fed by the 
proto-galactic disc, the BH accretes from the envelope. This struc- 
ture is dubbed a quasistar, since it resembles a (scaled-up) red gi- 
ant in structure, though it is p owered by accretion into a central BH 
JBegelman et al.ll2006l , l200a. hereafter BRA08). 



BRA08 found that for each black hole mass A/bh , there is a 
minimum envelope mass M, below which no hydrostatic solutions 
for the whole envelope can be obtained. This line corresponds to 
modest M,/Mbh ratios of a few tens. Above this line, the BH 
accretes at the Eddington rate for the whole mass (Mt = Mbh + 
M«). The BH can accrete at a super-Eddington rate (for the BH 
mass), because the energy released by the accretion is transported 
outwards through convection (and not by radiative transfer). After 
a few Myr, the BH reaches a maximum mass of ~ W^Mq, which 
provides a luminosity that unbinds the surrounding envelope. This 
physical picture has been confirmed by numerical simulations (Ball 
et al. 2011, hereafter BTZEll). Accretion will then proceed at a 
more "modest" rate through the proto-galactic disc. Accretion and 
mergers with other holes allow these "BH seeds" to evolve over 
cosmic times into SMBHs. 

An important characteristic of quasistars is that they are 
loosely bound, since they are radiation pressure dominated. A lumi- 
nosity that is close to the Eddington limit therefore has the potential 
to drive strong radiation-driven outflows. The presence of an out- 
flow, in turn, would modify the quasistar's internal structure, the 
BH accretion rate and the final mass of the BH seed that can be 
produced. Quantitatively, however, the extent and consequences of 
such winds are currently not known. 

In this paper, we aim to address the above questions. We calcu- 
late approximate solutions for a hydrostatic massive envelope pow- 
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ered by accretion, allowing for mass loss from its surface. We find 
that the presence of a wind allows for the existence of hydrostatic 
solutions for the envelope, even if the system is super-Eddington 
(for the total mass). With the wind, there is a competition between 
accretion and evaporation, and the BH can be effectively fed and in- 
crease its mass only when the envelope evaporation is sufficiently 
weak. This happens for very massive envelopes > IO^Mq, whose 
wind reaches a maximum mass loss, driven by all the available lu- 
minosity. 

We find that under the above conditions, a BH seed of > 10'' 
solar masses can be assembled in less than 10* yr. In particular, to 
obtain BH seeds of lO" - 10^ A/q, an envelope of mass > IO'^A^q 
is required. We will show that, as a consequence, the present sce- 
nario predicts the presence of massive BH seeds at high redshift 
{z > 10) only in halos of least pa a few 10® Mq. 

The paper is organized as follows. In i|2l we describe our 
model and carry out analytical estimates. In ij3] we describe the 
numerical solution for the structure of super-Eddington quasistars 
with winds, and we present our results in S|4] We discuss the impli- 
cations of our findings to the modeling of supermassive BH forma- 
tion in 



2 THE "QUASISTAR" STRUCTURE 

The system we consider is heuristically illustarated in Fig. [T] 
Crudely, it comprises of five primary regions. From the inside out 
they are: 

(i) A core, where an embryo BH is fed through a convectively 
dominated accretion disc. 

(ii) A convective envelope (which contains most of the quasistar 
mass). 

(iii) A porous atmosphere, where convection is inefficient. 

(iv) An optically thick wind. 

(v) A geometrically thin (at least on these scales) proto-galactic 
disk. 

The first four regions comprise the quasistar, the structure of which 
we solve here. 

The last region is that of the disk which feeds the quasitar. For 
reasonable galactic accretion rates, less than a few tens of Mq yr^ ^ 
the disc is locally sub-Eddington at the radius of the quasistao 
Since the disc remains slim, we can assume that mass exchange 
between the disc and the quasistar happens on a narrow equato- 
rial region. Moreover, the accreted mass is probably redistributed 
within the quasistar over a hydrostatic timescale, which is, as we 
will comment later, much shorter than both the accretion and the 
evaporation timescales. Thus, we can treat the disc just as a mass 
source term for the quasistar mass budget. Only in one case con- 
sidered here, for the extreme rate of 300 A/q yr^^, can the disc be 
geometrically thick and the interaction between the accretion and 
the wind from the quasistar surface may be dynamically important. 
Therefore, the results obtained in this case have to be regarded as 
indicative only. 

In the following sections, we describe in more detail the differ- 
ent components and construct an analytical model that will allow us 
to delineate the region in the A/* — A/bh parameter space, where 



^ Using the results of §2.2 for the radius r, and mass Af* of a n = 3 
polytropic stai; one can show that for M < 50 Mq yr~^, the accretion 
will release less than Lsdcb that is, M < LEdd/(G'M,/2r,). 



it is possible to grow massive BHs. These BHs can then serve as 
seeds for the SMBHs observed at the centre of galaxies. First, how- 
ever, we will state and justify the primary assumptions we employ 
to construct our model. 

For our purposes, it is not necessary to model the hydrody- 
namics of the innermost accretion flow, which by radial extension 
and mass is a negligible fraction (~ 10~^) of the quasisitar. The 
theoretical luminosity expected from this "hypercritical" regime of 
accretion is used as a source term in the energy equation for the 
envelope. Nevertheless, we do consistently consider the accretion 
rate which in turn depends on the envelope structure, as we explain 
in the next section. 

It has been shown that numerical modeling of quasistars is 
sensitive to the choice of the inner envelope radius (BTZEII). 
However, in absence of a physically sound alternative, the most 
natural choice remains the Bondi radius of the inner accretion disc. 
We therefore adopt it, as it facilitates the analytical calculations and 
the comparison of our results with previous numerical studies. 

It is possible to show analytically and numerically that a hy- 
drostatic envelope exists only if it is much more massive than the 
BH (BRA08, BTZEII). Therefore, we seek solutions satisfying 
this condition. On the other hand, there is an upper limit on the 
mass of quasistars, given by stability considerations. Such radiation 
dominated structures (with an almost null t otal binding en ergy) are 
in fact subject to the pulsation instability. iFowlerl ( 119661) showed 
that even a small degree of rotation can stabilize stars with mass 
< 10* A/q. Although the envelope is mainly pressure supported, 
quasistars are indeed expected to be (slowly) rotating since they are 
fed by pre-galactic accretion discs. We therefore assume A/* ~ 10* 
as an upper limit for the envelope mass. 

Rotation would also slightly squash the envelope at the poles. 
Nevertheless, it is not trivial to predict the exact rotational pro- 
file expected inside the quasistar. As we shall see, the quasistar 
will be almost entirely convective. We can therefore expect an az- 
imuthally dependent rotation rate, as exhibited in the sun. How- 
ever, this behavior is still not fully understood. As a consequence, 
we will model the quasistar while neglecting the effects of rotation, 
that is, while assuming spherically symmetry, as in previous works. 

2.1 The core 

In the innermost region (hereafter "the core", in analogy to stars), 
the embryo BH is fed through an accretion disc which trans- 
fers mass from the surrounding hydrostatic and much more mas- 
sive envelope. The extent of the core is the sphere of gravita- 
tional influence of the BH, the radius of which is the Bondi ra- 
dius Tb = GMbh/2c^, where the local adiabatic sound speed is 
Cc = Cs(rb) — \/4Pc/3pc. Here pc and Pc are respectively the 
matter density (p) and total pressure (P) at that radius. 

The small angular momentum in the system is such that 
very close to the last stable orbit, the flow is rotationally sup- 
ported and viscous processes allow the release of a luminosity of 
Lbh ~ TjMBnc^, with a radiation efficiency of 77 = t?o.i x 0.1. At 
larger radii, the disc has a "thick" geometry JNaravan & Yilll994l : 
lAbramowicz et alj|l995l : iBIandford & BegelmaJ|l999h. which is 
radiatively inefficient due to its high density JBegelmanI 1 1 9791 : 
iBegehnan & Meiej[T982l : lBlondinlll986l) . Since the radial convec- 
tion can be shown to be important, the innermost flow is therefore 
modeled as a convection dominated disc (CD AF, Stone et al. I999|; 
llgumenshchev & Abramowiczlll999l : lOuataert & Gruzinovii2000i") . 

At steady state, all the forces in these discs scale with grav- 
ity. Thus, Cs ex r~^'^ and a "modified hydrostatic" balance holds. 
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Figure 1. A heuristic description of tlie system we solve. See detailed description in ^ 



Moreover, contrary to classical advection dominated discs, CDAFs 
are characterized by a net (and constant) ou tvyard flux of en- 
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ergy, which is possible on ly if p c 
lOuataert & Gruzinovll2000h . 

These scalings imply that the maximum luminosity that con- 
vection can transport within the flow. 



A I i 

47rr pCg 



(1) 



is constant with radius. This limit for the carrying capacity of con- 
vection is set by the requirement that the convective cell motion 
should remain subsonic. That is, the convective velocity can be at 
most the local sound speed Cg, otherwise shocks will form and dis- 
sipate the motion over distances much shorter than the pressure 
scale height. Therefore, the energy released at small radii can be 
transported by convection up to ~ rb, if Lbh SJ I/c,max- We thus 
assume 



aic.n 



(2) 



with a $; 1. The parameter a not only accounts for convective 
transport inefficiency but also for possible leakage of energy (e.g., 
by jets close to the hole) which would decrease the supply of energy 
into the envelope. 

Since I/c,max is constant in the flow, it can be indeed evaluated 
anywhere not too close to the BH. However, for reasons that will 
be clear in the next section, it is convenient to express Lc,max at rb. 



iBH = a\/^7r(GAfBH)'p?^'r- 



(3) 



where we took into consideration that this hot flow is radiation pres- 
sure dominated, Pc — aT^ /'i, where a is the radiation constant. 



2.2 The hydrostatic region 

In a quasistar, most of the mass and the radial extent are occupied 
by a highly convective envelope. The convection efficiently trans- 
ports the accretion energy flux to large radii (~ 100 AU) where 
it powers a wind, after having crossed a thin hydrostatic radiative 
layer. In the following estimates, we will assume that the mass and 



the radius are constant in this radiative layer and equal to those of 
the entire convective envelope (M, and r, respectively). 

The structure of a strongly convective hydrostatic envelope 
can be described by an n = 3 polytrope (while employing the clas- 
sic Lane-Emden solution). These massive stars have a small fixed 
ratio between the gas and radiation pressure of 

-1/2 



/3fa7x 10~''(M,/10''Mq)" 



< 1. 



(4) 



Since the central density and pressure profiles are rather flat with 
radius, we can use the central polytropic values to estimate the den- 
sity and temperature at the Bondi radius, and thus evaluate Lbh (as 
given by eq.[3}. 

Scaling the envelope mass as Ah = m, A/q and the central 
temperature as Tc —Ta 10® K, the envelope radius is given by 



5.8 X 



inl2 1/2 

10 mj 



T." 



(5) 



while the central density is 
= 1.3 X 10" 



l/2rTn3 



Pc = i.a X iU 'm^-'-Ig gem (6) 

(e.g., see lHovle & Fowleilll963l). 

Furthermore. IJoss et alj ( 119731) have shown that wherever the 
density is high enough, the local luminosity can never exceed the 
local Eddington limit Ledd.r ~ 4:TTcGM{r)/K(r), as calculated 
with the enclosed mass M{r), and the local opacity K(r). This 
is because convection will be excited, and it will advect a large 
enough flux to keep the system sub-Eddington. However, there is 
an upper limit to the luminosity that can be transported by convec- 
tion at each location (eq.[T). Outside the Bondi radius, the density 
and temperature are nearly constant such that Lc,max oc r^. There- 
fore, the convective carrying capacity increases with radius near 
the centre. At larger radii however, both p and T decrease steeply, 
forcing a progressively higher fraction of the luminosity to be trans- 
ported by diffusive radiative transfer. Eventually Lbh = ic,max- 
This location, which we call the "transition radius" (rti), marks the 
base of the radiative layer. We search for solutions for which this 
region is in hydrostatic equilibrium, or in other words, solutions in 
which a radiative atmosphere is present. 

As atmospheres approach Ledd,r a plethora of instabilities 
can set in, whether the atmospheres are only Thomson scattering 
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20 lai). wh ether they ha ve magnetic fields ( lAronslll992l : 
iGartimie 1998; Be gelmanI 2002h. or even complex opacity laws 



dGlatzell 



19941 : IPapaloizou et al.lll997 i. which are probably not rel- 



evant for the primordial gas in quasistars). Once the instabilities 
grow to become nonlinear, they create low density channels for 
the photons to esca pe, and effecti vely suppress the radiative force 
exerted on the gas IShavivlll998B . Such a reduction in the opac- 
ity is the only known way to ex plain th e super-Eddington nature 
of classical novae and ?7-Carinae dShavivlllOOOa . 2001b) . Although 
we are far from having a complete understanding of the nonlin- 
ear state, we ca n account for it with a reduced opacity Kcff (e.g., 
iDotan & Shavi v 201 1). The functional form of this ejfective opac- 
ity HicS will be further elaborated in i]3.3l The Eddington ratio is 
defined as F = I//I/cdd,r(«:), where L(r) is the radiative luminos- 
ity. In the radiative atmosphere L{r) = Lbh. 

With a reduced opacity it becomes possible to have a static 
radiative layer, if we further impose the approximate condition that 
the luminosity at the transition radius equals the ejfective Eddington 
limit, then 



1.4 X 10^* Kt/ m, erg s" 



(7) 



where Str is the effective opacity at the transition radius (fi:(rtr) = 
Ktr), normalized to the electron scattering value, k^s ~ 0.35 cm^ 

In the radiative atmosphere, the luminosity satisfies Lbh ^ 
Ledd,r(Keff), Up to the radlus where the gas clumps become op- 
tically thin such that they cannot funnel the radiation anymore 
into the rarified regions. As the clumps become transparent, the 
opacity therefore recovers its microscopic value (kcH — >■ k) and 
Lcdd,r(«:cff ) — >■ icdd,r(K), implying that above the radiative atmo- 
sphere, there is a radius beyond which no hydrostatic solution can 
exist. This causes the quasistar to lose mass through a wind. In the 
whole section, we will assume that the wind sonic point is close 
to r,, i.e. Ts = r, (see section [J!4l for a detailed discussion of the 
wind solutions). 

Following BRA08, we combine eq.[3]and eq.|6l and by equat- 
ing the two expressions for Lbh (eqs.[3]and[7](, we derive an ex- 
pression for the central temperature: 



rn 1-7 in4 / - ^-2/5 -4/5 7/10 t^ 

Tc ~ 1.7 X 10 (a/s;tr ) "'■bh "^* K. 

If we insert it into eq.[5]and eq.[6l we finally get 



^ o c inl4 I ~ \2/5 4/5 

r, « 3.5 X 10 (oKtr) mgjj 
and 
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(8) 



(9) 



(10) 



Pe ^ 6.4 X 10-^° (a^tO"'^' m-^n'^mT^ 
respectively. 



2.3 The "no-solution line" 

We now proceed to calculate the line, in the Mbh — A/, param- 
eter space, below which no hydrostatic envelope is possible. We 
term this boundary as the the no-solution line. As we shall see, it 
corresponds to low envelope masses (M*/Mbh ~ 10). 

Without mass loss, solutions with progressively smaller total 
envelope masses (for a given BH mass) have progressively lower 
photospheric temperatures (BRA08). We find numerically (section 
3) that the presence of a wind does not alter this conclusion. As the 
envelope mass decreases, the photospheric temperature will eventu- 
ally reach the recombination temperature, below which the opacity 



drops precipitously. Any further decrease in the envelope mass re- 
duces drastically the photospheric opacity, causing the photosphere 
to reside further inside the wind, where densities are higher. This 
continues until the photosphere reaches the base of the wind, that 
is, Tph — Tb- When the envelope mass is small enough to al- 
low this to occur, the photospheric temperature will be given by 
LBH/(47rrJ) 



<^T^Y, 



, which can be evaluated to be 



Tph = 1.1 X 10^ C"° "~'/' ^bT mV^'' 



K, 



(11) 



where a is the Stefan-Boltzmann constant. 

When Tph — Ts, the presence of an entirely optically- thin 
wind does not affect the structure of the hydrostatic envelope and 
we should recover the solutions derived for a purely hydrostatic 
quasistar. Under that assumption, BRA08 have shown that — in 
analogy to the Hayashi track for red giants and protostars — no hy- 
drostatic solution can be found below a minimum photospheric 
temperature Tmin, which is just a factor of 2 lower than the re- 
combination temperature. It was shown that Tmin varies at most 
by a factor of two around 4500 K for BHs in the range between 
A/bh = 1 to 10'' A/q (see their fig. 3). 

Assuming a constant floor temperature and solving for Tph = 
Lmin, we find that the "no-solution line" is given by 

20/7 



11 ir,4 4/7 -9/7 / Li 
7n» = 1.1 X 10 Q ' K, ' ' 



V4500K 



mBH\*/'' 



100 



(12) 



We note here that, while BRA08 link the existence of this limit 
to a steep drop in opacity in the radiative layer, BTZEl 1 argue that 
it is a general feature of polytropic models and not one linked to 
the adopted opacity law. Both papers however, agree that this limit 
exists and that it corresponds to AL/A/bh ~ 10, which is what we 
obtain with eq.[T2](see Fig. 2). 

2.4 The wind and the evaporation strip 

It should be stressed that even if a hydrostatic envelope can exist, it 
does not guarantee that a quasistar can "live" long enough for the 
BH to grow inside it. In the presence of strong winds, we should 
also require that the accretion time-scale, 

A/bh ALbh 



^BH — —■ — -^ 7 ^, 

Mbh Lbh/??c^ 
is shorter than the evaporation time-scale. 



to 



Af, GMi 



Mw e Lbht-* 



(13) 



(14) 



Here, we normalized the mass loss rate to the maximal possi- 
ble rate, obtained when all the accretion luminosity is used to 
drive the wind out of the gravitational potential well, A/„ 

iBH/(«o'=c/2) « LBHr./GAL., 



^w.max 



AL„ 



14 



o 2/5 ~-3/5 /mBH\4/^ /m, \-l/5 _^ 



(15) 



where Uosc is the escape velocity. That is, we have written that 

2 

eLbh, (16) 



1) 



with an efficiency factor e = Af„/AL„,max =S! 1. 

Using the expression for the quasistar radius (eq. |9]l and rj 
O.lryo.i, weget 



, 11/5 -9/5 

— Ki 4.2 X 10"^ '"* "^BH 

iBH ' e(Q«:tr)^/^»7o.i ' 



(17) 
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The above ratio, at the no-solution hne (eq. I12t . is smaller than 
unity for Mbh < IO^Mq, even if we assume that only 1 percent 
of the luminosity has been used to accelerate the wind (i.e., e = 
10-2). 

Let us now fix A/bh and explore what happens as we increase 
the envelope mass. First, we note that the strongest dependence is 
on the envelope mass. Therefore, as we increase A/*, fcv/teH in- 
creases and it will eventually be greater than unity. Also, the radius 
of the star r, decreases and the depth of the potential well from 
which matter should escape increases: v'^^c oc M^, j k^^ . 

Since the mean density M,/r^ increases as well, it is rea- 
sonable to assume that, as the envelope mass increases, the radial 
distance between the sonic radius and the photospheric radius in- 
creases. This implies that an increasing fraction of the luminosity 
is used to pull the gas out of the potential well (i.e., e — >■ 1). At that 

point, Afw ^ A/w.max. 

We anticipate here one of our numerical result (elaborated in 
@, that the condition icv/ieH > 1 is found above the e = 1 limit 
for all BH masses we consider (A^bh < IO^A/q). Therefore, the 
equal evaporation/accretion line is given by 



10^(afi:ti 



x2/ll 2/11 /mBHN^/" 

' "^"-^ \ 100 



(18) L^ = 



There is a relatively wide strip in the A/* — A^bh parame- 
ter space between the no-solution line (eq.ll2b and the equal-time- 
scale line (eq.llSt, that we dub the evaporation strip. Only above 
this strip is it possible "to grow" black holes inside quasistars. 
Therefore, we will refer to the tcv/iaH ~ 1 line as the threshold- 
growth line, and to the region above as the growth region. 



2.5 The quasistar Evolution 

Let us now assume that the envelope is constantly fed at a rate of 
a few Az/q yr^^ from the pregalactic disc. From eq.[T6] we deduce 
that, in fact, the mass of the envelope will likely not increase, since 
the mass loss will (over)compensate the feeding. An upper limit for 
the BH mass that can be reached is therefore obtained by setting 
the envelope mass to be constant in time. The quasistar will then 
move rightwards in the A/, — Mbh place and it will eventually 
cross the threshold growth line. Afterwards, the envelope mass will 
decrease abruptly until the quasistar will hit the no-solution line 
and the envelope will disperse. 

The maximum mass for the BH is therefore given by eq. [TS] 
This implies that it is necessary to have Af, > 10® A-Zq in order to 
form a BH seed, more massive than ~ 10"^ Mq . We will discuss the 
implications of our findings in ^ after having verified this picture 
with numerical calculations. 



such that the threshold-growth line lies in the e — 1 region. How- 
ever, in order to determine the actual extent of the evaporation strip, 
it necessary to calculate numerically the structure of the quasistar. 



3.1 The innermost accretion region 

As was the case in the analytic analysis, we do not model the ac- 
cretion flow. We assume that a luminosity Lbh given by eq. |2]is 
transported through the quasistar. 



3.2 The Hydrostatic Convective Envelope 

Above the accretion zone resides the hydrostatic region of the en- 
velope. It is composed of two parts. The inner part is the convective 
zone, in which the energy emitted by the accretion is convected out- 
wards. Convection is an efficient process of energy transfer only as 
long as the luminosity Lbh is lower than the maximal convective 
luminosity, Lc,max, given by eq\T\ 

The hydrostatic envelope is described by the equation of hy- 
drostatic equilibrium, 

ldP_ _GM{r) 
p dr 



(19) 



where M{r) is the total mass contained within radius r. 



M{r) = A/bh + / 4Tvr pdr, 

the equation of state, 

k,Tp ^ 1 ^^,^ 



P^P^ + Pr^ 



/imp 



(20) 



(21) 



where fcb is the Boltzmann constant and /imp is the mean mass 
per particle. Given our primordial composition and that in most 
of the quasistar T > W* K, we adopt fi — 0.5, independently 
of radius. The last assumption does not affect our global results 
for two reasons. First, quasistars are radiation pressure dominated. 
Second, the pressure is important only over the hydrostatic region 
and the base of the wind, where the temperature is always above 
the ionization temperature. The temperature gradient is given by 
the adiabatic gradient. 



dT _ r2 - 1 T dP 

Hr ^ r2 P d7' 

where r2 is the adiabatic index: 

_ 32 - 24/3 - 3^2 
^ ^ 24 - 18/3 - 3/32 ' 

and P = Pg/P. 



(22) 



(23) 



3 NUMERICAL SOLUTION 

The full modeling of a quasistar, including the radiatively driven 
wind, cannot be carried out entirely analytically. In particular, one 
carmot derive the wind strength and its optical depth as a function 
of M, and A/bh. In part, this is because the stellar "boundary" 
conditions at the sonic point depend on the characteristics of the 
wind, and vice versa. 

In the previous section, this ignorance was buried in the un- 
known behaviour of the parameter e, and of r(rtr) that enters the 
equations through kh. We simply assumed that their evolution was 



3.3 The Hydrostatic Radiative Region 

According to common wisdom, objects with a mass M* cannot 
shine beyond their classical Eddington limit, I/Edd ~ 47rGAf«c/K, 
since no hydrostatic solution exists. In other words, if objects do 
pass I/Edd, they are highly dynamic. They have no steady state, 
and a huge mass loss should occur since their atmospheres are then 
gravitationally unbound. Thus, persistent astrophysical objects can 
pass I/Edd but only for a short duration corresponding to the time 
it takes them to dynamically stabilize, once super-Eddington con- 
ditions arise. 

For example, this can be seen in detailed ID numerical simula- 
tions of thermonuclear runaways in classical nova eruptions, which 
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can achieve super-Eddin ston luminositie s but only for several dy- 
namical timescales (e.g.. Istarrfieldl 19891) . However, once they do 
stabilize, they are expected and indeed do reach in the simulations, 
a sub-Eddington state. Namely, we naively expect to find no steady 
state super-Eddington atmospheres. This, however, is not the case 
in nature, where nova eruptions are clearly super-Eddington for du- 
rations wh ich are orders o f magnitude longer then their dynamical 
timescale IShaviv"2001b). This is exemplified with another clear 
super-Eddington object — the great eruption of r^-Carinae, which 
was a few time s above its Eddington luminosity for over 20 years 
( IShavivll2000bl) . 

The existence of a super-Eddington state can be naturally ex- 
plained, once we consider the following: 

(i) Atmospheres become unstable as they approach the Edding- 
ton limit. In addition to instabilities that operate under various 
specia l cond it ions (e.g. , Phot o n bubbles in strong magnetic fields, 
AronsI 1 19921 : lOammiel 1 19981 : iBlaes & Socr"atesll200ll : Igegelman' 
2002 , or s-mode instability under special opacity laws, Glatzel 
19941 : IPapaloizou et alj|l99'7[). two instabil ities operate in Thom- 






son scattering atmospheres JShavivluOOlaf) . It implies that all at- 
mospheres will become unstable already before reaching the Ed- 
dington limit. 

(ii) The effective opacity for calculating the radiative force on 
an inhomogeneous atmosphere is not necessarily the microscopic 
opacity. Instead, it is given by 



off 



(Fkv) 

{F)y 



(24) 



where ( )y denotes volume averaging and F is the flux JShavivl 
The situation is very similar to the Rosseland vs. Force 



opacity means used in non-gray atmospheres, where the inhomo- 
geneities are in frequency space as opposed to real space. For 
the special case of Thomson scattering relevant here, the effective 
opacity is always reduced. 

For r ^ Fcrit, the effecti ve opacity can be desc ribed by an 
empirical parametrization (e.g.. lDotan & Shavivll201 ih 



«^<=ff^'^( i-f^) T' 



(25) 



where Fcrit is the critical F above which inhomogeneities are ex- 
cited. F rom theoretical considerations, we take Fcrit ~ 0.8 dShavivl 
l2001a J). In our work we choose B = 1, which implies A — 0.16. 

Finally, we assume that the relation between the effective Ed- 
dington factor Fefi = -£//-£/Gdd,r(«;eff) and the classical Eddington 
factor r = L/Lcdd.r('i) is empirically given by 



A 



r for r < Fcrit 



for F > Fcrit, 



Fefl = 

Fcff = F for F < Fcrit. (26) 

As was found in lDotan & Shavivl ( 1201 ID . the solutions we find in 
21 are relatively insensitive to the exact choice of A and B. 

The central temperature of a quasistar is of the order of 10® 
K (see eq. 8) and it decreases slowly in the convective envelope. 
Therefore, the temperature in the hydrostatic part of the quasistar is 
high enough (^ 10* K), that the opacity is mainly due to electron 
scattering, and therefore approximately constant, k ~ Kcs . 

Unlike the hydrostatic part, the temperature in the wind can 
drop below 10'* K, especially near the photosphere. This means that 
hydrogen recombines and opacities due to "bound-free" transitions 
become important. Since we are modeling objects that are formed 



ou t of pristine or mildly po lluted gas, we use the opacity computed 
bv lMaver & Duschll ( 120050 . In particular, we use an analytical fit 



l + (T/8 X lO^K) 



nscm'g ' 



(27) 



(BRA08), which captures the rapid drop in opacity for T < 10*K. 
The temperature gradient is given by the radiative one with the 
above effective opacity. 



dT_ 
dr 



IQ-nacr'^T^ ' 



(28) 



We note that in the atmosphere, the radiative luminosity is con- 
stant and L — Lbh, because work has not been done on the gas. 
In fact, the atmospheres of super-Eddington objects, effectively re- 
main sub-Eddington while being classically super-Eddington. This 
is true only as long as the inhomogeneties comprising them are 
optically thick. This condition will break at some point where the 
density is low enough and the wind will ensue. 



3.4 The Wind 

When inhomogeneities in the atmosphere becomes optically thin, 
the opacity returns to its microscopic value, and a net outward force 
acts on the gas. We will show in the following, that the nature of 
the outflow and the mass loss rate ultimately depends on the ratio 
of gas to radiation pressure and therefore on the quasistar mass Af* 
(eq.©. 



3. 4. 1 Two-fluid solution 

In a wind, the radiation energy can be transported by diffusion and 
by advection, so that the total energy rate is L + Afw4Pr/p, where 
the latter is the advected enthalpy, while L is the diffusion luminos- 
ity, 

LW = -?i^^^^. (29) 



dr 



Sup 



In a low mass quasistar — we will quantify it shortly — the 
mean density is too low for advection to be efficient and the main 
mechanism for radiation energy transfer in the wind is diffusion. 
The radiation can thus be described as an additional external force 
acting on the gas, besides gravity. For this reason we refer to this 
solution as a "two-fluid" wind. 

Using the continuity and Euler equations in steady state, it is 
possible to show such trans-sonic wind solutions will be charac- 
terized by having an isothermal sonic point which coincides with 
the critical point. The "isothermal" sonic point is the radius where 
V — Cyj = ^Pg/p (with Cw being the isothermal speed of sound), 
while the critical point is the radius where the net forces balance 
each other. 

Using the above definition of the sonic point, the mass loss 
rate can we written as 



Af„ = 47rrsPwCw = const. 



(30) 



where, p^ is the density at the sonic point. 

In our case, the critical point will be located where the net 
radiative force on a gas element balances gravity, namely, where 
Fcff ~ 1. This will happen when the inhomogeneities comprising 
the porous atmosphere become optically thin such that they cannot 
reduce the opacity anymore. This point will also mark the outer 
extent of the hydrostatic part of the quasistar, rs ~ r« . 
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Based on the fact that the inhomogeneities are the resuh of ra- 
diative hydrodynamic instabilities which operate on length scales 
comparable to the density scale height in the atmosphere, it is pos- 
sible t o estimate the average density at the sonic point JShavivl 
UOOlbh . Using this density, the mass loss can be estimated to be 



M„ 



WLbh 



(1 - r.-) 



(31) 



where Fs — r(rs) and W is a dimensionless wind function JShavivl 
uOOlbh . In principle, W can be calculated ab initio only after the 
nonlinear state of the inhomogeneities is understood. This however 
is still lacking as it requires elaborate 3D numerical simulations 
of the nonlinear steady state. Nevertheless, one can compare the 
above prediction to observed mass loss rates in classical novae and 
in giant eruptions of L uminous Blue Variables, and find that W ~ 
5 — 15 (Shaviv 2001b). We shall take a value of 10. By comparing 
eq. [16] and eq. [3T] we have an expression for the efficiency factor 
e = W^«esc/(2ccw) (1 - rr') ^ 5uL/(cCw). 

As mentioned before, when the mass of the quasistar in- 
creases, the compactness of the star increases as well. As a con- 
sequence, an increasingly high efficiency e is necessary for the pre- 
dicted mass loss to be pushed to r ^- cxd by the available lumi- 
nosity. When the required efficiency satisfies e > 1, which means 
Vesc ^ \Jc^c/y\!, a time independent wind solution is then not 
possible, because the wind would have to stagna te at a finite radius. 
The r egime is then that of "photon-tired winds" dOwocki & GavlevI 
[T997h . 

The time-dependent behavior of photon tir ed winds was nu- 
merically investigated bv lvan Marie et alj ( 120091) . It was found that 
shocks form between infalling material and the outflowing wind. 
This forms a layer of shocks in which there is a large kinetic flux, 
but without the associated mass flux. In particular, they find that the 
mass loss from the top of the layer of shocks is reduced to e < 1: 
ef»min[0.2r^*',0.9]. 

Therefore, in general, the wind mass loss is given by 



M„ 



M„ 



o.2^^■^o.9, 



M„ 



Mw 



(32) 



Beyond the sonic point, acceleration to, say, v ~ 2cw happens 
in a geometrically thin layer of a few scale heights. In the following 
supersonic regime, the thermal gradient is no longer important with 
respect to the external forces. From this point on the gas velocity is 
described by 



dv 
dr 



GMt [1 - F(r)] 



(33) 



The last equation that is needed to describe the wind is that of 
energy conservation, which for r > Vs reads 



Lbh = L{r) + M„ 



1 2 



GMt GMt 



= const, 



(34) 



where the advection term has been consistently dropped. 

We are now in the position to quantify the regime in which this 
"two-fluid solution" can apply. The requirement that the gas flows 
through a sonic point, beyond which F > 1 is equivalent to the 
requirement that the advective term must not reduce the luminosity 
below the Eddington value by r^. 



Mw < Lbh 

P 



Le 



(35) 



This criterion will be slightly different depending on the mass loss 
regime, which determines A/w, but in both cases it can be ex- 



pressed in terms of a minimum gas to radiation pressure ratio. In 
non-photon-tired winds, (eg.lBlt, it becomes 



/3>4W— . 
c 



(36) 



We remind here that Cw is the isothermal sound speed calculated at 
the sonic point. In the photon-tired regime, the criterion for having 
an isothermal critical point is instead. 



P> 



8F 
F-1 



^csc 



(37) 



Because generally one has both that Cw <^ c and Cw <^ foac, 
the two-fluid solution with an "isothermal" sonic point breaks down 
only for very low gas pressure fractions, which in quasistars implies 
a large envelope mass. Numerically, the line above which there is 
no isothermal critical point is denoted in fig. [2] One can see that ex- 
cept for the least massive black holes, this line resides in the photon 
tired regime. This is because the non-photon-tired criterion (eq.|36b 
is generally more stringent than the photon tired criterion (eq. l37l l. 
Eq. [36] implies a mass range 



M, 



IQf^Mp 



<0.3 







lO'^cm/s 



-1/2 



(38) 



3.4.2 Adiabatic Winds 

At low enough gas pressures, no isothermal critical point can de- 
velop since the luminosity in the fluid frame of reference is sub- 
Eddington. Instead, the radiative energy is mostly advected with the 
flow. One can then consider the radiation and gas to flow together 
as one fluid, with an effective adiabatic index of the combined fluid 
(which is very close to 4/3). 

In this regime, the quasistar will develop a slow outflow above 
the porous atmosphere, enough to keep the flux sub-Eddington, but 
still sub-sonic. The flow would only become supersonic further out, 
where the escape velocity is similar to the local speed of sound, 
resul ting in an adi abatic wind (i.e., the classical wind solution for 
stars. lParkeJiigeoh . 

To see that such a subsonically outflowing layer must be 
present, let us assume for a moment that it is not. In such a case, 
an adiabatic "single fluid" wind will develop with a sonic point lo- 
cated where the polytropic part of the quasistar ends (eq. [9]l. To 
find the energy required to drive a wind from that location, we 
need to consider that the polytropic object has P — Kp*''^ with 
K = (M/87r)^/^7rG, that its radius is given by eq.[9] and that for 
a thermal wind, v ~ Ucsc at the sonic point. The ratio between this 
energy and the energy released by the black hole is then 



Lwind 2-Krlpv^vl^^ 



12 






> 1. 



(39) 



Clearly, the required energy is much larger than the amount avail- 
able. This implies that the thermal wind cannot be accelerated from 
where the polytropic part of the star ends. Therefore, there must 
be a layer with an outflow which is fast enough to keep the lumi- 
nosity sub-Eddington (relative to the gas), but slow enough to be 
subsonic. Only further out, where the escape speed is significantly 
slower, will the flow become supersonic. 

Since the amount of energy left in the flow at the sonic point 
is much smaller than the energy it had deep below, it is clear that 
almost all the energy is used to pull the material out of the gravita- 
tional potential well, and very little will be left as kinetic energy in 
the thermal wind, or in the form of observable luminosity. 
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Thus, because the system is hmited by the gravitational well, 
adiabatic winds are characterized by the similar mass loss rate 
Afw,max of photon tired winds, albeit in a configuration which is 
markably different from before. 



3.5 Numerical Methods 

Practically, we perform the numerical calculation of a quasistar as 
follows. 

(i) Choose a black hole mass, A^bh- 

(ii) Choose an inner total pressure, Pc. 

(iii) Guess the gas to total pressure ratio at the inner radius, /3c = 

With the above values, one can calculate the inner temperature, 
density, adiabatic speed of sound, Bondi radius and luminosity by 
respectively using 



3x103 



Tc = 



3(1 - /?c)Pc 



-,1/4 



Pc = PcP, 



kTc ' 



Cc = \/4-Pc/3pc 



rb 



GMbh 
2c2 • 



(40) 



(41) 



(42) 



(43) 



The luminosity Lbh is given by eq.|2]with a = 1. In fact, Lbh is 
evaluated at 5 x r-b, to avoid the density cusp generated by the point 
mass potential of the BH (see eg. 119b . 

Using the above values, we integrate the equations of stellar 
structure for a convective envelope (eqs. [T9] - l22t until Lbh = 
I/c.max. When Lbh > Lc.max, we use the same equations except 
that the temperature gradient is instead given by eq. [28] This inte- 
gration is carried out up to Keff ~ 1- This location, where FcS ~ 1, 
is the wind sonic radius rs in the two-fluid regime (sec. 13.4. Ft . From 
this point outward, we integrate the wind equations (eqs. [33] [32l 
and 134b. with the following initial conditions. At r = rs we take 
the value of the gas isothermal sound speed c(rs) and tempera- 
ture r(rs) given by the hydrostatic solution and we set the ini- 
tial wind sound speed, velocity and temperature to be Cw ~ c{rs), 
Vw = Cw and T„ = T{rs). Assuming mass continuity at r^ with 
the wind rate given by eq.[32l we derive the initial wind density 
Pw ~ Af„/47rrsCw. The way we connect the hydrostatic solution 
to the wind uses the fact that the acceleration happens in a very nar- 
row radial region, and we approximate it to happen at one radius, 
r = rs. 

Finally, the wind equations are integrated up to the photo- 
sphere r — Tph, where the optical depth is r = 2/3. At rph, the 
temperature r(r-ph) should be such that one recovers consistent 
surface conditions for the radiation field: 



r(rph) 



47rr2^c 



1/4 



(44) 



To enforce it, the value of /3c (our only free parameter) is iterated 
and, for each successive guess, the quasistar structure is recalcu- 
lated until eq.|44]is satisfied. 

Since the adiabatic wind conditions are satisfied above the 
photon tiring limit (Fig. |2}, we do not need to calculate specifi- 
cally that solution, since the mass loss rate is going to be similar to 
the one of the corresponding photon tired solution. 
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Figure 2. The chai'acteristic regions in the M» — A/bh pai'ameter space. 
Separated by black lines are: (1) The region with no hydrostatic solution 
(bottom) (2) The evaporation strip, where the envelope evaporation time- 
scale is shorter than the black hole growth time-scale (3) The growth region, 
where evaporation is less important. The short-dashed black line represents 
the thi'eshold quasistar mass above which the winds are photon-tired. The 
blue lines represent lines of equal Mm [Mq yr~^]. The dotted lines are 
the horizontal equi-BH accretion rates, in units of AIq yr"'^. The long- 
dashed red lines show the analytical predictions for the threshold-growth 
line and the no-solution line. The analytical threshold-growth line assumes 
a constant F = 10, while the analytical no-solution line assumes constant 
r = 3, T„,in = 4500 K and the numerical K^g. 



4 NUMERICAL RESULTS 

In this section, we first discuss the different regions in the M, — 
A/bh parameter space for the structure of a quasistar Then, we will 
follow its temporal evolution towards the formation of massive BH 
seeds. 



4.1 The "no-solution" region 

Using the numerical analysis, we first confirm the existence of a 
region in the AI, — A/bh parameter space where no hydrostatic 
self-gravitaing envelopes can be found. The region corresponds to 
AL,/ALbh ^ 18, as depicted in Fig|2] In this region, the hydrostatic 
layer between r^, where the supersonic accretion begins, and rg, 
where a the supersonic wind starts, is simply too thin geometrically 
to be stable (with rs/r^ ^ 2). 

The boundary of this "forbidden" region in the A/* — A/bh 
parameter space corresponds to envelopes with a photospheric tem- 
perature of 4500 < Tmin < 6000 K, in agreement with BRA08. 
In particular, we find that Tmin increases from 4500 K to 6000 K 
between A/bh = 5OA/0 and A/bh = IO'^A/q; then it decreases 
down to 5400 K at 10^ Mq- For a given A/bh, Tmin is the min- 
imum photopheric temperature with which a quasistar can shine, 
which corresponds to the minimum possible envelope mass, ac- 
cording to Tph oc ALJ/^" (eq.lHJ. 

We compare our numerical results with the analytic prediction 
given by eq.[T2](lower long-dashed red line). It assumes a constant 
Fs = 3 and Tmin ~ 4500 K. The good agreement reinforces the 
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Figure 3. Wind mass loss and Eddington ratios as a function of tlie envelope 
mass, for a fixed BH mass of A/bh = IOOMq. The Eddington ratios are 
calculated at r^ (Fg) and at the photosphere r-ph (Fobs = r(rpij)). 



validity of our numerical model. Despite the more sophisticated 
modelling of the atmosphere which allows for mass loss, we re- 
cover the expected behaviour in the limit of low envelope masses. 
As the envelope mass decreases, for a given BH mass, r-ph -^ r-g 
and the wind becomes completely optically thin and dynamically 
unimportant. 

4.2 The Evaporation Strip 

At the no-solution line, quasistars are already very windy despite 
the relatively modest Fs, with a loss of a few to ten thousand A/q 
yr^^, as shown by the lines of constant Afw in Fig. |2] (solid blue 
lines). There, the vigorous mass ejection causes the evaporation 
time-scale to be shorter than the accretion time-scale, tcv < ten- 

The mass loss intensifies as a given BH increases its enve- 
lope mass, Afw oc M^'"^ (Fig. [3j. It will eventually reach the 
photon tiring limit, after which it increases at a much slower rate 
Mw oc M~'^^^r^^^ oc M°''^ (see eq.[T6]and Fig.[3j. This can also 
be seen in Fig. |2l where above the short-dashed line, the lines of 
constant Afw become almost vertical. Correspondently, the evapo- 
ration time-scale first decreases as icv oc Ai^"'^, then, once in the 
photon tiring regime, it increases as tcv oc A/^'®. 

The accretion rate onto the BH, instead, always increases lin- 
early with the envelope mass. Therefore, the accretion time-scale 
goes as fBH oc M~^. Because of this faster and monotonic evolu- 
tion, it eventually equals the evaporation time-scale inside the pho- 
ton tiring region. 

The existence of an "evaporation strip", where the envelope is 
lost to the wind rather than feeding the BH is thus confirmed by 
our numerical calculations (Fig. |2]l. The extent of the strip, given 
by the line of equal time-scales (tcv = ian) is well in agreement 
with our analytical prediction (eq.[T6l upper red long-dashed line), 
where we use a constant Fs — 10. As we explained in section [23] 
we call this line "the threshold growth line", since it indicates the 
minimal envelope mass within which a given BH can significantly 
increase its mass by accretion before the envelope is lost. 



4.3 The evolution of a Quaslstar 

So far, we focused on single equilibrium solutions. However, the 
structure depicted in Fig. 1 is in constant evolution. On the one 
hand, the envelope loses mass to the wind and to the BH. On the 
other, it gains mass, accreting from the proto-galactic disc. Since 
the dynamical time-scale \/ R^/GM, ~ 20 m-^^m, yr, is 
always much shorter than both the accretion and the evaporation 
time-scales, we can construct the quasistar evolution as a succes- 
sion of equilibrium states. 

For this description to be valid, we require that the thermal 
timescale will be shorter than the evolutionary timescale. To esti- 
mate it, we should note that quasistars are to a good approximation 
n = 3 polytropes, in which the internal energy almost compensates 
the gravitational binding energy, leaving a small contribution from 
the finite gas pressure. Namely, 



Eu 



-m 



: l3k2GM^'^ p^^ , 



(45) 



where ^2 = 0.639. Since the thermal timescale is roughly tth 
£tot/iEdd, we find using eqs. 4, 7 & 10 that 



t\h < Etot/L-Edd ~ 2m 



'^''^Btyr- 



(46) 



As we shall see, this timescale is shorter or at most comparable 
to the evolutionary timescale. 

Once we are allowed to consider steady states, we can envi- 
sion two scenarios. In the first scenario, we assume a large amount 
of mass collapses to form the quasistar, after which it equilibrates 
and evolves while the proto-galactic disk continues to be accreted. 
In the second scenario, we assume that the quasistar forms with a 
relatively low envelope mass, and it grows while the accretion rate 
is larger than the wind mass loss rate. 

As study cases for the first scenario, we follow the fate of a 
BH of IQOMq embedded within envelopes of initial mass of 10®, 
10^ and 7 x 10^ Mq (Fig.|4j. In each case, the quasistar is born in 
the accretion zone. We considered several different accretion rates 
onto the envelope, Afacc = 1, 3, 10, SOA/q yr~^, and we found 
very similar tracks, since the change in Af« is mainly due to wind 
losses, for the range of envelope and BH masses its track passes 
through. We therefore plot one examples with Afacc = 10 Mq 
yr-\ 

As the BH mass increases by accretion, the quasistar moves 
towards the right of the plot, until it encounters the evaporation 
strip after a typical time ~ 10^ yr has elapsed (see the temporal 
evolution in Fig.[5]l. This evolution in the accretion zone is so short, 
that the envelope loses only about 25% of its initial mass. However, 
as soon as the quasistar enters the evaporation strip, the envelope 
is stripped off at such a high rate that the quasistar plunges down 
vertically in the diagram, reaching the no-solution line in a further 



lO'' 



10'' yr. 



The main result here is indeed that the presence of the evap- 
oration strip limits the configurations in which an embryo BH can 
grow to masses > 10'' A/q. In particular, it is necessary that a qua- 
sistar "is born" above the line of threshold growth limit, with an 
envelope of at least l(f Mq. The initial envelope mass determines 
the final BH mass. Empirically, for an initial BH mass of 100 Mq 
and A/facc = 10 Mq yr~^, we find 



^'^-■' - ^6° (life 



Mp 



(47) 



where AfeH.f is the final BH mass and Af,^ is the initial enevelope 
mass. This relation holds for any reasonable accretion rate, Afacc < 
lOOA^Q yr"\ and any initial BH mass A/bh,! > IOOA^q. For a 



10 Dotan et al. 



much lower BH mass, the evaporation rate may be lower than the 
accretion rate and there may be a dependence on the initial BH 
mass, likewise, for much higher accretion rates. 

In the second scenario, we assume that the quasistar is formed 
through a more gradual accretion process. For the "low accretion 
rates", at least several dozen times the Eddington rate, a bare BH 
can handle the accr eted mass with a disk l ike solution, having op- 
tically thick winds (Dotan & S havivll201lb . However, the BH will 
not be able to handle higher accretion rates with this configuration. 
This will necessarily puff up the system to form a "minimal" qua- 
sistar. From this point, we wish to integrate its evolution. We expect 
the envelope mass to increase until the wind will compensate for the 
accreted matter. That is, it will evolve by increasing the mass of the 
cocooned BH, until the system will reach the no-solution line. 

Two such tracks are plotted in Fig. [4] and the temporal evolu- 
tion is given in Fig.|5] These quasistars are formed in the evapora- 
tion strip with i\fBH,i = 25AIq. Here we find that the final mass 
of the BH is determined by the maximal envelope mass which in 
turn depends on the mass accretion rate. 

For Mbh.i = 25 Mq, M,,! = 1.4 x IO^Mq, we calculate 
numerically relations between the accretion rates and the final BH 
mass. For accretion rates < 20A/q yr^^, which do not allow the 
quasistar to leave the evaporation strip, we get the relation 
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(48) 



For higher accretion rates which place the Quasistar above the 
growth line, we get 



A^BH.i 



550 



Ma, 
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(49) 



,lOM0/yr, 

Only in this last case, the final BH can be larger than 10^ Mq, but 
only for Mace ^ 30 Mq yr~^. However, it is a weak function 
of Mace, and already extremely high accretion rates (> 100 Mq 
yr~^) are needed to counterbalance the envelope evaporation and 
grow A/bh with at least a few 10'^ Mq. 

In addition, we note that the case of an intial BH with a 
few tens of Mq is a favourable configuration in which a quasis- 
tar formed in the evaporation strip can somewhat grow its BH. As 
we move towards the right within the evaporation strip, the winds 
become more and more vigorous (see Fig. 2) and the envelope is 
stripped off progressively more rapidly, leading eventually to con- 
ditions that will not allow any BH growth. 

We thus conclude that, for plausible galactic accretion rates 
(~ 1 — 10 Mq yr~^), quasistars should be necessary formed above 
the evaporation strip with an envelope of at least IO^'Mq, in order 
togetMBH > 10^ A'/q. 

4.4 Observational appearance 

During the growth of the BH, quasistars are very dim objects, since 
most of the radiation energy is converted into kinetic energy for the 
wind. The photons that do manage to escape from the photosphere, 
make up a modest luminosity which is smaller than the Eddington 
(see Fobs in Fig. [U. 

When the quasistar enters the evaporation strip, however, the 
luminosity starts increasing fast as the envelope is being blown off. 
Eventually, Fobs ^ 1 beyond the photon tiring limit, but this Super- 
Eddington phase is very brief (< 10'^ yr), where most of the time 
is in fact spent at the photon tiring limit. 

Another interesting point is the fact that because of the wind. 
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Figure 4. Quasistar evolutionary tracks. The blue lines describe the evo- 
lution of quasistai's which accrete from the pre-galactic disc at a rate of 
10 Mq yr~^. The upper three ones, have an initial BH mass of 100A/q, 
and initial envelope masses of 10®, lO'^ and 5 X lO''^ Mq, respectively. 
The quasistar with 10® A^q spends 1.5 X 10"^ yr in the accretion zone and 
a comparable time, 6 X 10'' yr, in the evaporation strip. The aforemen- 
tioned time-scales for Af* = IO'^A^'q (and M, = 6 X lO'' Mq) are 
1.3 X 10^ yr (600 yr) and 1.5 x 10^ yr (500 yr) respectively (see Fig.|5). 
In all cases, the mass evolution of the envelope is governed by the wind 
loses. In the BH growth zone, the quasistai's lose only Si 25% of their ini- 
tial mass, while all of it it is lost in the evaporation strip. The final BH 
masses are: 750AfQ, 1.3 X IO'^A^q and IO^A/q. The lower blue line is 
for a quasistar fomed in the evaporation strip with Mbii,i = 25 Mq and 
Af,.i = 5.7 X 10^ Mq, while the red-line is for the same intial conditions 
but with Afacc = 300 Mqyi'^ . In these last cases, the final BH mass is 
determined by accretion rates. There temporal evolution is shown in Fig.|5] 



the objects never have a high effective temperature. As a conse- 
quence, any cosmological redshift will imply that the objects can 
only be observed in the infrared band. Therefore, we can conclude 
that although they were very brief transients, they can in principle 
be detected by the forthcoming James Webb Space Telescope. 



5 DISCUSSION AND CONCLUSION 

In the present work, we considered the state and evolution of qua- 
sistars, where a massive envelope enshrouds a black hole. We found 
that an important characteristic of these objects is the continuum 
driven winds which they accelerate. This gives rise to an evapora- 
tion strip in the envelope mass - BH mass parameter space, where 
quasistars evaporate. BH growth takes place only for large enough 
Af, for which the wind mass loss is limited by photon tired winds. 
We then considered two evolutionary scenarios. In the first 
scenario, the quasistar is initially formed above the evaporation 
strip, by first having a large gas cloud collapse, then form an in- 
ternal black hole, and subsequently evolve. In the second scenario, 
a seed BH is first formed, and then it accretes gas at large rates. We 
found that only the first scenario can generate large seed BHs. This 
is because once a BH is formed, unrealistically high mass accre- 
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Figure 5. The temporal evolution of the BH (solid lines) and envelope 
(dashed lines) masses for three quasistars. The magenta lines describe the 
case with Mbh,! = IOOMq and Af»_i = lO'^ A/q. Here the quasis- 
tar is formed in the growth region (see Fig. |4) and reaches AfBH,f = 
1.3 X 1O*M0, in 2.8 x 10^ yr. The blue and black lines describe the 
evolution of quasistars with an initial BH mass of 25Mq and M»^i = 
5.7 X W^Mq, which is accreting at a constant rate of lOM© yr~^ and 
300 Mq yr~ ^ , respectively. The final BH formed in this scenario heavily 
depends on the accretion rate. It has a mass of 44OA/0 in the first case and 
almost 5.7 X 10^ Mq in the second case. Both quasistars live for f» 1 
Myr Note that the thermal timescale in the latter case is comparable to the 
evolutionary timescale, and therefore it should be considered cautiously. 



tion rates are required to evolve the quasistar out of the evaporation 
strip. 

These results put strong constraints on the dark matter ha- 
los in which massive SMBH seeds can form. Our first scenario, 
on which we no w concentrate, follows the picture put forth by 
iBegelmanI yOlOj), whereby the progenitors of Quasistars are su- 
per massive stars that formed as the consequence of the high infall 
rate of hydrogen-cooled gas at the centre of dark matter halos. The 
life-time of these stars is set by the thermonuclear timescale for 
burning their hydrogen core. At the Eddington limit, it is '--^ 2 Myr, 
independent of mass. After that, the core may collapse into a BH. 
Once the BH starts accreting from the envelope, the feedback from 
the released luminosity sets up the structure we have investigated 
in this paper. 

Our results constrain the initial mass of the quasistar, and thus 
of the supermassive star, that can grow massive SMBH seeds. A 
quasistar with an initial mass ~ 10^ Mq (or « IO^Mq) can form 
a BH of ~ IO'^A'/q (or ~ lO^AfQ respectively, as can be see in 
Fig. 4). This initial mass should be accumulated in less than ~ 2 
Myr, which requires accretion rates greater than a few AIq yr~^. 

Dark matter halos accrete matter through their virial radius at 
a rate of 



A/i 



4.6 {z + l)?o' 



f- 



M. 



M 



©yr 



(50) 



JNeistein et alj2006n . where Afh is the mass of the dark matter halo 
that hosts the quasistar and (z + l)io — {z + 1)/10. 



S ince the timescal e for dark matter and gas accretion is the 
same JDekel et al.ll2009h . 



Mace = /b X Mb 



(51) 



it is also an estimate for the accretion rate that feeds the quasistar, 
under the assumption that almost all gas can be funnelled towards 
the halo centre. In eq. | 5 1l /b ~ 0.17 is the cosmological baryon 
fraction JKomatsiJIlOllh . 

In fact, it is not clear if such high accretion rates do not un- 
dergo fragment ation at parsec scales and if so, how much gas turns 
into stars (e.g.. lShlosman & Begelmanll 19891 : lLeviiJl2007l) . This is 
known as the "active galactic nuclei (AGN) fueling problem", since 
we observe them shining with a luminosity that implies ^ IMq 
yr~^ . However if, like it must happen in AGNs, t his problem is not 
severe jWise et al.ll2008l : lBegelman & Shlosmanli2009 ) such flows 
can reach the proto-galactic centre. 

Under this assumption, we can thus connect the accretion rate 
needed in order to form a quasistar of a certain mass and at a certain 
redshift, with the mass of the host halo through eq.|5T] However, we 
should also assume an efficiency factor, since not all gas /b x Mb 
in the halo can be used to form a quasistar. Therefore, the minimum 
halo mass that can host a given quasistar with mass A/, is 



A/h 



M 







7x 10^ f M« 



{l + z)l,^ VlO«A% 



,6 X 10** 



M. 



106 A/q 



,(52) 



where we assumed in the second term that no more than 1% of the 
total amount of gas is used for the formation of the quasistar. The 
first term comes directly from eq. [sT] 

Eq.[52]implies that at z = 10, quasistars with M* > 10® Mq 
can be found only in dark matter halos with A/h > 10^ Mq. In 
particular, for z > 10 "massive" SMBH seeds of A/bh > lO** Mq 
need host halos with Mb > 6 x 10^ Mq. 

The very bright quasars observed at z ~ 6 have masses of typ- 
ically 10® A/q, estimated using the Eddington argument. Moreover, 
seeds of 10* — 10^ A/q can grow by accretion at the Eddington 
rate in ~ 0.5 Gyr. Therefore, such SMBH seeds may have formed 
as late as z ~ 10. At this redshift, the comoving number density 
of ha los with Mb ~ 10^^A/q (calulated with the Sheth & Torm3 
1 19991 formalism) matches the observ ed comoving nu mber density 
of those bright quasars, 10"®Mpc"^ jFan et alj200lh . Those halos 
are indeed above the threshold set by eg. 1521 

Of course, halos formed at diiferent redshifts can lead to the 
observed high redshift quasars and to the observed local galaxies 
with SMBHs. To properly derive the consequences of eq. [52] one 
should thus follow a "merging tree" evolution of halos and also of 
the hosted black holes. This is beyond the scope of this paper, but 
it will be addresses in a forthcoming work. 

Finally, our results show that, unfortunately, Quasistars will 
not be easy to observe. They are relative rare objects which shine 
at super-Eddington rates only for a short period, prior to their final 
evaporation. 
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